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Abstract 

We have completed the next-to-leading order perturbative QCD corrections 
to the virtual-photon exclusive differential cross sections for heavy quark pro- 
duction in deep-inelastic electron-proton scattering, i.e. e + P^Q + Q + X. 
Using these results, we have computed distributions which are sensitive to cor- 
relations among the heavy quark, the heavy antiquark, and the associated jet. 
Some predictions for charm and bottom heavy quark production at the electron- 
proton collider HERA are presented. 
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1 Introduction 



Order a s QCD corrections to structure functions containing heavy quarks and to 
single heavy quark inclusive distributions in deep- inelastic electroproduction (i.e., 
7* (<?) + Pip) ~> Qipi) + X where X stands for any final hadronic state allowed 
by quantum-number conservation and Pip) is a proton of momentum p ) were 
recently published in jl[] and respectively. By combining the next-to-leading 
order (NLO) heavy quark structure functions with the corresponding light-quark 
structure functions ||, the heavy quark content of the nucleon has been studied as 
a function of Q 2 = —q 2 and x = Q 2 /2p • g Q. Event rates for charm production 
integrated over bins in x and Q 2 relevant to HERA data have been calculated in || . 

To further the study of deep-inelastic electroproduction of heavy quarks we have 
recalculated the virtual-photon-parton cross sections of pj in an exclusive fash- 
ion. This enables us to study the single and double differential distributions and 
correlations among all outgoing particles in the virtual-photon induced reaction 
"Y*+P^>Q + Q + X with X = or 1 jet and to easily incorporate experimental 
cuts. A short letter, containing a study of invariant mass distributions of heavy- 
quark-heavy-antiquark pairs, has been published |J. In this paper we present all the 
calculational details and additional distributions. The transverse and longitudinal 
photon components are treated separately and the latest CTEQ3 parton densities 
0], consistent with the newly released HERA data B, are used in the kinematic 
regime appropriate for production of charm and bottom quarks at HERA. We make 
our predictions at fixed values of Q 2 (> 8.5(GeV/c) 2 ) and x (> 4.2 x 10~ 4 ). 

Our treatment of electroproduction at HERA extends the existing studies of in- 
clusive QCD corrections in the virtual-photon channel |l]], inclusive QCD corrections 
in the real-photon channel M, and exclusive QCD corrections in the real-photon 
channel fTofl , allowing for an extensive comparison with present and future experi- 
mental data. Heavy quark electroproduction has already played an important role 
in the determination of the gluon distribution function in the proton at low x. In 
fact, a study of the photoproduction and the Weizsaker- Williams electroproduction 
of charmed mesons has just appeared [11 1. Production cross sections of charm and 
bottom quarks are also relevant in the proposed study of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements through the rare decays of D— and B— mesons 
and the DD and BB mixing [12]. Our study of deep-inelastic electroproduction is 
motivated by our desire to avoid complications associated with the parton densities 
of the photon which must be incorporated when Q 2 is small. Therefore our results 
should give a cleaner test of NLO perturbative QCD. By examining distributions 
and correlations that are trivial at lowest order (for example, the azimuthal angle 
between the heavy quark and the heavy antiquark) one directly tests the hard scat- 
tering cross section that is predicted by perturbative QCD and factorization. Here 
we should note that this is a fixed order perturbative calculation and suffers from 
the same problems as all NLO calculations. Hence there are regions in phase space 
where it will break down. For example, in the above mentioned azimuthal angle 
distributions, if one looks at the prediction too near the back-to-back configuration 
there will be an extra enhancement of the cross section due to multiple soft gluon 
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emission which is not included in our NLO result. 

We stress that here we only consider extrinsic heavy quark production, involving 
Bethe-Heitler and Compton production from massless partons. For a discussion of 
intrinsic production, where the heavy quark is considered to be part of the proton's 
wavefunction, see Brodsky et al. fl3|| . A variable flavor scheme which joins the 
extrinsic heavy flavor production picture at /U p h ys <C m with a light mass parton 
density picture at ^ p h ys 3> m has been discussed by Aivazis et al. [jT|]. By comparing 
the fixed flavor scheme calculation of |jj with the variable flavor scheme of [|14|] it 
is concluded in that the former yields the most stable and reliable results for 
F2(x, Q 2 ,m 2 ) in the threshold region ( i.e., Q 2 < 10m 2 where m is the mass of the 
heavy quark). 

In our exclusive computation we use the subtraction method which is based 
on the replacement of divergent (collinear or soft) terms in the squared matrix el- 
ements by generalized plus distributions. This allows us to isolate the soft and 
collinear poles within the framework of dimensional regularization without calculat- 
ing all the phase space integrals in a spacetime dimension n 7^ 4 as usually required 
in a traditional single particle inclusive computation. The method has been used 
in the context of electron-positron annihilation [|l(|, hadroproduction of Z° boson 
pairs [17 1, hadroproduction of jets [18], hadroproduction of heavy quarks [19[| , pho- 
toproduction of heavy quarks flO| , and hadroproduction of a W boson plus a photon 
The essence of the method is nicely described and compared to the also popu- 



lar phase-space slicing method [21] in the introduction of the paper by Kunszt and 



Soper [18|. The paper of Mangano et al. |H| contains many useful details. In this 
method the expressions for the squared matrix elements in the collinear limit ap- 
pear in a factorized form, where poles in n — 4 multiply splitting functions and lower 
order squared matrix elements. The cancellation of collinear singularities is then 
performed using the factorization theorem p^| . The expressions for the squared 
matrix elements in the soft limit appear in a factorized form where poles in n — 4 
multiply lower order squared matrix elements. The cancellation of soft singularities 
takes place upon adding the contributions from the renormalized virtual diagrams. 
Since the final result is in four-dimensional space time, we can compute all relevant 
phase space integrations using standard Monte Carlo integration techniques [^3] and 
produce histograms for exclusive, semi-inclusive, or inclusive quantities related to 
any of the outgoing particles. We can also incorporate any reasonable set of exper- 
imental cuts. Our computer code has no small phase space slice parameters and 
the parameters defining the generalized plus distributions may be tuned to give fast 
numerical convergence, which is an advantage of using this subtraction method. 

We therefore have a new, and more general, calculation of the scale indepen- 
dent (Wilson) coefficient functions (or scale independent partonic cross sections) 

^hi^Jli^) 1 Cfci^'O) an d ^ii^iC) defined in Q, as functions of r\ = s/Am 2 — 1, 
where s is the square of the c. m. energy in the virtual-photon-parton system, and 
£ = Q 2 /m 2 . We checked the n and £ dependence of these scale independent coef- 
ficient functions and the x and Q 2 dependence of the hadronic structure functions 
F2(x,Q 2 ,m 2 ) and Fl(x,Q 2 ,m 2 ), obtained after convolution of the coefficient func- 
tions with the light quark and gluon densities in the proton, against the results in 
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We comment on this comparison below. The results for the single quark in- 
clusive transverse momentum and rapidity distributions, dFk(x,Q 2 ,m 2 ,pt)/dpt and 
dFk(x, Q 2 , m 2 , y)/dy (k = G,L) respectively, were also checked against the results 
in ||] to make sure that all our computer programs yield consistent results. 

In sect. 2 we introduce our notation used in the calculation of the cross sections 
for the off-shell photoproduction of heavy quarks. We closely follow the notations 
of |IJ and [ 10 1 . In sect. 3 we discuss the analytical results needed for the calculation 
of the numerical results. The latter are presented in sect. 4. We concentrate on 
those plots which are dependent on the emission of an extra parton (or jet) and 
only show single differential distributions. Our conclusions are presented in sect. 
5. Kinematics and other useful results are relegated to appendix A. Appendix B 
contains a discussion of the soft limit of the matrix elements. Appendix C contains 
a discussion of the collinear limit of the matrix elements. We have included these 
appendices to help with the identification of the various pieces of our computer code. 



2 Notation and Born cross section 



In this section we introduce our notation and define the quantities which we will 
calculate. We then review the Born level results. The kinematics and other useful 
results are given in appendix A. Following we denote the parton subprocesses 
contributing to heavy quark production in deep-inelastic electron-hadron scattering 
by (see fig. 1) 

7*(g) + ax(kx) -» Q( Pl ) + Q(p 2 ) + a 2 (k 2 ) + ■■■ + o i (fe / ), (2.1) 

where the photon momentum is space-like ( q 2 = —Q 2 < ), ai(ki) stand for 
massless partons (k 2 = 0), and Q (Q) is the heavy (anti) quark (p 2 = p\ = m 2 ). We 
consider the partonic structure tensor 

W, u = + / , 2 #a 17 [£M,(2)M:(2)dT 2 +£M M (3)ilC(3)dr 3 + •••], 

(2.2) 

where each term in the square brackets corresponds to a perturbative expansion in 
a s . The terms in the structure tensor are the flux factor l/2s', the initial degree 
of freedom average 1/2(1 + S gai e/2) where S gai is the Kronecker delta function, the 
color average K ai ^, the 2 to j body amplitude M(j), the j body phase space dFj, 
and the sum over initial and final degrees of freedom S. The color average factor has 
the specific values K gi = l/(N 2 — 1) and K qi = K gi = 1/N, where the number of 
colors is N = 3 for SU(3). We work in the context of dimensional regularization with 
space time dimension n = 4+e. Since we integrate over the azimuthal angle between 
the plane containing the incoming and outgoing leptons and the plane defined by 
the incoming parton and outgoing heavy quark the partonic tensor ( |2.2| ) only has 
two terms which can be written as 

W, u = dffr{-g^ + ^>j (2.3) 
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V J ^ ^2 _9l/ J I J ^ T L > ' 

Using the projection operators g^ v and fc^i the partonic cross sections can be 
written as follows: 

*"> ~ -Itttt^-- {2A) 

d<TL = -^jfcffcfrvV (2.5) 
V s ) 

We work with da G and dc^ for convenience and later use these to get the transverse 
partonic cross section dox via the relation 

1 1 

da T = da G + 2 \ + e /2 dcTL ' 

For each of the three partonic channels a\ = g,q,q we expand the corresponding 
parton cross sections in terms of the number of outgoing particles in reaction ( |2.1[) : 

da hai =da% i+ da% i+ --- (2.7) 

where i = G,L and 

d^l=C hai Mt\j)dV r (2.8) 

Using eq. (|2.2| ) we find the partonic cross sections are related to the squared matrix 
elements by 

M lU) = "(IwEWi^Kt^i-j), (2-9) 

with 

c -'» = Qtt^) °< (e) (^srnb/s*^) • (2 ' 10> 

where aa(e) = 1 and ax,(e) = 2(1 + e/2). Note that 

a, g = 2C i ,(l + e/2)C l - g , (2.11) 

with the color factor Cp = (N 2 — 1)/2N = 4/3 for SU(3). The expressions for 
MqQ) and M^(J) were evaluated using the symbolic algebra program FORM 



We begin our analysis by considering the lowest order process contributing to 
reaction fl2,l|) which is shown in fig. 2. Namely, 



Y(q) + g(h)^Q(pi) + Q(P2). (2.12) 

The projections of the matrix element for this reaction are calculated according to 
eq. (2.9). We calculate in n = 4 + e space-time dimensions. Both the O(e) and 0(e 2 ) 
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terms are needed for a correct treatment of renormalization and mass factorization. 
Summing over initial and final spins 



Mf (2) = 8g 2 e 2 e 2 H NC F B kjQED , 



(2.13) 



with k = G,L and 
Bg,qed 



2 c' 



ui h 2q z s 4m z s 
-T + — + T— + -, 



1 I 

m s 



2^2 



+ 



2m q 



i\2 



+ e- 



1 (s') 2 *V q 4 
-1 + ^^ + — ±- + 



m 2 q 2 (s') 2 



Bl,qed 



-V 

'^2 



tllil tiUi t\Ui 

m 2 (s') 2 ' 



Atiui 



(2.14) 
(2.15) 



The kinematic variables are defined in appendix A. Here g and e are the strong and 
electromagnetic coupling constants respectively and en is the charge of the heavy 
quark in units of e. The lowest order cross section now follows directly from eq. 



3 Order a s corrections 

In this section we describe the calculation of the NLO corrections to the lowest order 
gluon fusion cross sections which were derived in the last section. As mentioned in 
the introduction, we proceed via the subtraction method wherein one extracts the 
soft and/or collinear singularities without integration over final state particles Q 
and a2 as done in the previous single particle inclusive calculation |l]], ||. This 
allows us to plot correlations between final state particles which will be presented 
in sec. 4. We first discuss the gluon-bremsstrahlung reaction in detail and then, in 
the following subsection, the light quark contributions. The kinematics of the three 
body final states are discussed in appendix A, which also contains the definitions of 
the variables used here. 

3.1 7*# Channel 

The gluon-bremsstrahlung reaction, shown in fig. 3, is 

7*(g) + g(h) -» g(h) + Q( Pl ) + Q( P2 ). (3.1) 

We begin by writing the partonic cross sections as 

dtr® = C i>g M?(3)dT 3 

= C i , g f?(x,y,6 1 ,9 2 )dr 3 , (3.2) 

where 

fi(x,y,9 1 ,9 2 ) = t\u' - q 2 s 5 /s)M[(3), (3.3) 
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dT 3 

t> («' - q 2 s 5 /s) ' 



(3.4) 



We suppress all arguments of ff other than x,y, 61,82 for notational convenience. 
These are the variables we integrate over to calculate the partonic cross sections for 
reaction ([O]). The subtraction method proceeds by multiplication of the squared 
matrix element, and division of the phase space, by invariants (or combinations 
thereof) that vanish in the soft and/or collinear limit. The invariants are chosen 
to cancel the singularities in the squared matrix elements due to the propagator (s), 
thus giving finite functions ff(x, y, 9%, 82) which can safely be implemented in a 
computer program. The modified phase space dT 3 is expanded using generalized 
plus distributions allowing extraction of soft and/or collinear poles in n — 4 which 
are cancelled as usual. 

In this particular the photon is virtual, and the heavy quarks are mas- 

sive, the only collinear singularity comes from the collinear emission of a gluon from 
the incoming gluon as shown in fig. 3(d). The squared matrix element Mf(2>) has a 
1/(1 + y) singularity when y — > —1, therefore it suffices to multiply by t' oc (1 + y)(l - 
to get a result that is finite as y — ► — 1. Similarly, the emission of a soft gluon cor- 
responds to a 1/(1 — x) 2 singularity in Mf (3) as x — > 1 so we choose to multiply by 
an additional factor of (u 1 — q 2 s^/s) oc (1 — x) to get a result that is finite as x — > 1. 

We now proceed to show the cancellation of the singularities and derive the 
renormalized partonic cross sections. Substituting the expressions for if and u' 
found in appendix A and using eq. ( |A.12 ) we find 

-l+e/2 



-HNdT 



(5), 



-l+e/2 



7T 



: I - x)- 1+t (l - y z )- l+t ^dysm t 6 2 d62 , 

(3.5) 

where all quantities are defined in appendix A. We now begin to extract the poles, 
by expanding the (1 — x)~ l+e piece in terms of the generalized plus distributions 
defined in appendix A. Using eq. ( |A.l §| ) one finds 

-l+e/2 



da. 



(3) 

i,9 



-C itg HNdTf( 

7T 



-l+e/2 



J\-l+e/2 



5(1 -x) 



1 



+ 2 In + 2e ln z [3 



+ 



(1 
2 



dysm e 6 2 d92f!(x,y,9 1 ,62) 



C itg HNdvf( 



J\ -l+e/2 



-l+e/2 



vr \ s / V 1 

(1 - y 2 y 1+e/2 dy sin e 9 2 d6 2 f?(x, y, 6 U 8 2 ). 



+ e 



ln(l 



x 



(3.6) 



The 5(1— x) /e piece in the first term is a soft singularity now regulated dimensionally. 

x) = 5(1 — x)dT 2 dx while in the second piece we 
expand (1 — y 2 )~ l+e l 2 using eq. (|A.17| ) to find, 



In the first piece we use cZr 2 5 ^(l 



da. 



(3) 
'-.</ 



d(7 W + + da p 



«,9 



where 



^2 = -Ci, g HNdT 2 (s') 



-l+e/2 



7T 



*,3 



-l+e/2 



1 



+ 21n/3 + 2eln 2 /3 



(3.7) 



fi(0i), (3-8) 
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with 



and 



da) 



f!(Oi) = J 6(l-x)dx(l-y 2 )- 1+e/2 dysin e 9 2 de 2 fi(x,y,e 1 ,e 2 ), (3.9) 

/ln(l — x) 



>-<j 



-a„HNdr 2 5 \ s ')- 1+£/2 

7T 



/x -l+e/2 



1 — X 



+ e 



V l 



x 6(1 + y) [- + lln2oj) dysm e e 2 d9 2 ff(x,y,e 1 ,e 2 ), 



p. 
(3.10) 



and 



da) 



1 4) 



-C iig HNdTf{s')- 1+ ^ 2 (-^ 



7T 



1 — X 



+ e 



ln(l — x) 
1 — x 



l + y 



+ 



i-y 



dysm e 9 2 d9 2 fi(x,y, 9 1: 9 2 ). 



Pi 

(3.11) 



According to the factorization theorem 22] the cross section for scattering of a 
virtual-photon off a hadron H of momentum p is 



J 



(3.12) 



where the sum runs over all partons in H and we have suppressed the polarization 
index temporarily. Here (fti/ni^ A 4 /) is the probability density for finding parton i 
in hadron H with momentum fraction between ^ and £ + d£ at scale /ij. Noting 
that the infrared safe d&i is independent of the external hadron, we are free to set 
H = j, where j represents a parton. Then the partonic cross sections satisfy 



daj(ki) J dai(xk\)(f>i/j(x, fi 2 )dx. 



Up to first order in a s 



<t>i/j(x,Hf) = <%<5(1 - x) + — Pij(x)- + fij(x,n f ) 



(3.13) 



(3.14) 



where Pij(x) denote the Altarelli-Parisi splitting functions [25| and the functions 
fij(x,/j, 2 ) depend on the mass factorization scheme chosen. In the MS scheme 



fij(x,fj,j) = Pij(x) {je ~ hi47r + In fiff/J 2 



(3.15) 



where 7^ is the Euler-Mascheroni constant. In this paper we choose the mass 
factorization scale [if to be equal to the renormalization scale /x so that the last 
term in eq. ( 3. 15j ) is zero. Expanding daj and d&i in powers of a s gives 



daj(h) = daf{ki) + dof>{ki) + 



(3.16) 
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dffi(xfci) = da$ 0) (xki) + daf^x&i) + • • • . (3.17) 
Substituting eqs. (PI) , ( p5|) , (gig) , and ( gig) into eq. (gl3|) and equating powers 



of a s we obtain 



daf (fci) = daj 0) (fci), 



i 

where 2/e = 2/e + 7^ — ln47r. 

In the collinear limit 1/ -> -1 we show in appendix B that 



(3.18) 
(3.19) 



where 



and 



fi(x,-l,9 1 ,9 2 ) = ff(x,9 1 ) + f?(x,e 1 ,9 2 ), 
ff(x,9 1 ,9 2 )sm e e 2 d9 2 = 0, 



(3.20) 
(3.21) 



ff(x,0x) = 512vr 2 ^ 2e a 2 e 2 e 2 H NC A C F 



l-x 



x(l — x) + — 



s 

x 1 — x 

— + 

— X x 



B^QEDixh). (3.22) 



Performing the y and 6 2 integrations in eq. (3.10) using eqs. ( |3.20| ), (3.21), and 
flpl ) we find 



da. 



ln(l — x) 



(5) 



1 



1 e V 1 — x 
1 — x 



+ 2 



x(l — x) + 



1 — x 



+ 



+ 



p 

1 - X 



1 — X 



1 s ' 1 s ' w 
In — ~ + m h In — 



1 — x 



Bi,QED(xki). 



(3.23) 



Appealing to eq. ( |3.19| ) and reinstating the polarization index i we find 



d^+da^+da^ 



2^ 



| dx^(xA;i)P ra (x)| (3.24) 



where we have used dcj ^ = = 0. The Born level cross sections dafl are given 



i,q 



by eqs. and ( 2.13 ) and the virtual corrections denoted by da\* are discussed 
below. It is well known that EE 



2CU 



2C A 



(l-x). 

x 

(1 " ^ 



1 — X . 
H h x(l — xj 



+ 



x 

l-x 



X 



+ x(l — x) 



+ ( y Ca - |2>n, / + AC A In /3 ) S( 1 - ,• j . 



(3.25) 
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where Tp = 1/2, and nif is the number of light quarks. In the second line we 
have exchanged the usual plus distribution for the more general p distribution and 
defined $ = -y/1 — p. Cancelling the pole in eq. ( |3.23| ) by using eq. ( 3.24| ) gives the 
"factorized collinear" term 



da. 



i,9 



32C i;g a 2 s e 2 e 2 H NC A C F B i! Q ED {xk 1 )drf 
1 



1 



, s' , s' , u 
m — = + In h In — 

pr s 2 



+ 2 



(l-x) 2 + l + 

/ln(l - a 
V 1 -x 



(l-x)< 



(3.26) 



In the limit Q 2 = this result reduces to formula (2.45) of [[nj] when one identifies 



M§ } = 32Q gTra s e 2 e 2 ^NC fBg,qed- The remaining poles in eq. ( |3,24| ) cancel when 



the soft and virtual cross sections are added according to 



da lg 



A W l A ( s ) 
d(J i,g + d °i,g 



2 ni 



16Ci, g ajp- e e 2 e 2 H NC F dT 2 - —C A - -T F n lf + AC A \n(5) B^qed^ 



(3.27) 



The soft cross sections da\ s g \ which follow from eq. (|3.§|), are given in appendix B. 
The one loop order eg 3 diagrams for the virtual corrections to the reaction ( 2.1 2| ) are 
shown in fig. 4. The interference of these diagrams with the lowest order diagrams 
of fig. 2 were calculated in ]lj using the renormalization scheme of [26 and are 
available as FORM code. Using ( |3.271 ) we analytically checked the cancellation of the 
singularities. Also available from p]] as FORTRAN code is the analog of eq. (3.27) 
using the A prescription (27| for dividing the phase space into soft and collinear 
regions. The components of ( p. 27 ) in the A prescription and the subtraction method 
differ by finite terms. Therefore, we modified the original FORTRAN code by adding 
and subtracting appropriate finite pieces. In fact, the soft finite pieces to subtract 
off are given in eqs. (3.24) and (3.25) of |lj, and the pieces to add on are the finite 
pieces of eqs. ( B.16j ) and (B.17) of appendix B. Note that there is a typographical 

error in eq. (3.25) of Q; the sign of the In 2 r s term should be negative not positive. 

( f) 

Finally the three body contributions da\ are finite in n = 4 space-time dimen- 
sions and a short calculation from (|3.11| ) shows they may be written as 



da 



>-<j 



(l67T 2 ) 



Cinf3< 



1 



x ff(x, y, 9i,02)dxdy sin#id#id#2 



(3.28) 



where /?5 = i/l — 4m' 2 /s^. The finite functions /f(x, y, 0\, 62) follow from the def- 
inition (3.3), eq. (2.9), and the diagrams in fig. 3. Summarizing, we have the final 
result 

da i;g = dafl + da^ + da^ + d*% , (3.29) 

that is finite in n = 4 space-time dimensions and is to be used in eq. ( 3.12j ) to make 
predictions. Each of the last three terms in ( 3.29| ) individually depends on p and 
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However, the sum is independent of them as it must be because the decomposition 
( 3.29Q follows directly from eqs. ( A.17| ) and ( |A.18| ). This provides a strong check on 
our computer code. 



3.2 7*g and 7*5 Channels 

Analysis of the partonic reaction(s) 

7*(?) + QihMh)) - q(k 2 )(q(k 2 )) + Q(p x ) + Q(p 2 ) , (3.30) 

proceeds as above only this time there are no soft or virtual contributions making 
the analysis much simpler. This time we write the partonic cross sections as 

^2 = C iiq M9(3)dT 3 

= C l , q f!(x,y,9 1 ,e 2 )'dT;, (3.31) 

where 

f?(x,y,9 1 ,9 2 )=t'M?(3), (3.32) 

and 

dT^ = dT 3 /t'. (3.33) 
Replacing the divergent factor (1 + y)~ 1+e ^ 2 in ^3 by eq. ( [A. 16 ) one obtains the 



following decomposition: 

rh,. ,,, 

i,q i,q 1 i,q 



d° ( ?} = do\ c - ) + do\fl (3.34) 



with 



and 



d*t ] = -lc hq HNdv { i\ s 'Y /2 {^f 2 (1 - xY(i - y y/ 2 

x 5(1 + y) (- + Iulu) dy sin* 9 2 d9 2 f?(x,y,9 1 ,9 2 ), (3.35) 



dtzJJ = --C ls HNdTf{sr 12 f-V /2 (1 - xf{l - yY t2 



7T V S 



x (— J— ) dysm e 9 2 d9 2 f?(x,y,9 1 ,0 2 ). (3.36) 



As in the gluon-bremsstrahlung reaction (see appendix C) 

ff(x, -1, e u 6 2 ) = ff(x, 9 X ) + f?(x, 9 U 9 2 ) , (3.37) 

with 

f?(x,e 1 ,0 2 )sm € 9 2 de 2 = O. (3.38) 

Jo 
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In this case we find that 



Appealing to eq. (|3.19|) with 

P gq {x) = C F 



128TT 2 fi- 2t a 2 s e 2 e 2 H NC F (l + e/2) 



1 + (1 - x) 2 + ex 2 /2 



Bi,QED(xh 



(3.39) 



1 + 



one finds 



da 



(3.40) 
(3.41) 



with 



do - , 



»>9 



8C l ,,a2 e 2 e^iVC F S ii Q £Z) (xA;i)drf ) 
C l + (l-x) 2 



1 + 



(/) 



16vr 2 



s' s' a; 
In + In — + In - + 2 ln(l 
\i l s 2 



(3.42) 



C i , q p 5 f?(x,y,6 1 ,9 2 ) 



1 



1 + 2/ 



dxtiy sin#i<i#i(i#2- 



(3.43) 



The finite functions f?(x,y, 61,82) follow from the definition ( |3.32| ), eq. (|2 . 9[) , and 
the diagrams in fig. 5. In the limit Q 2 = 0, the "factorized collinear" cross sections 
for the quark channel daf "* reduces to formula (2.62) of [JO]] when one identifies 



g 7ra s e 2 e'j I NCFBG t QED and uses eq. fl2.11|) . This completes the analysis 
of the quark induced reaction. The analysis of the antiquark induced reaction is 
completely analogous. 

As the quark channel only contains collinear singularities, we will use it to illus- 
trate how the generalized plus distributions are implemented numerically and how 
the u> dependence disappears in the sum ( 0.41 ). To this end consider 



da, 



if) 



dyf{y) 



1 



1+y 



-l+LJ 



-1 



dyf(y) 



l + y 



+ 



/fo)-/(-D + 

-.1 i + y 



dyf(y) 

■l+w 

1 *w 

-x+u l + y 



l + y 



-1 1+9 



-\ i + y 



(3.44) 



where we have suppressed all indices and arguments of ff(x, y, 8\, 62) other than 
y. In the bottom line we see that the infinity encountered at the lower integration 
limit y = —1 is cancelled in the sum of the two integrals. In practice one can only 
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reasonably take the lower limit to be —1 + 5 where 5 ~ 10~ 7 in double precision 
FORTRAN before round off errors enter. None the less, the final results are stable 
with respect to the variation of 5 in the range 10~ 5 to 10~ 7 . The upper limit of 
the second integral gives a contribution /(— l)lno; which cancels against the lno; 

*■ (c— ) 

term in da\ . The first integral in the bottom line is commonly called the "event" 
and has a positive definite weight. The second integral plus the factorized collinear 
contribution is commonly called the "counter-event" and may have either positive 
or negative weight. The implementation for the gluon channel is similar but more 
complicated due to the presence of the soft-virtual terms. 



4 Results 

4.1 Virtual-photon-parton cross sections 

In this section we discuss the virtual-photon-parton cross sections and compare them 
with the results of M. The total cross sections are obtained by integration so that 



<7fc, 



da. 



k.i ■ 



(4.1) 



where daj-j are given as eqs. ( 3.29 ), and ( |3.4lD . Choosing a renormalization scheme 
where the heavy quarks in the gluon self-energy loops decouple in the limit of small 



momenta flowing into the loop [26], we express the perturbative expansion of the 
virtual photon-parton cross section in terms of scaling functions as follows: 



2 2\ 

7 ,m 



act* 



F(l) 



(77,£)ln— ^ 



(4.2) 



where 



1, (=^, 

4m 2 m 2 



(4.3) 



with s the square of the c. m. energy in the virtual-photon-parton system and /j, 2 

the mass factorization scale. Since the /^](77,C) depend on the electric charges of 
the heavy and/or light quarks we extract these charges and define new functions via 



<•(!), 



=2 Jl) 



.2 j(l) , 



(4.4) 



together with the corresponding formulae for the coefficient of the mass factorization 
term. Note that the o^ (r],£) vanish after integration to form the total photon- 
parton cross sections so they were not needed in [Q. However, now they must be 
retained when plotting distributions. Further, (rj, can only be evaluated at 
£ = provided the additional collinear divergence which arises when the photon 
goes on-mass-shell is mass factorized. As we only consider Q 2 > 8.5(GeV/c) 2 we do 
not perform this factorization although it was performed in |ffl] and checked against 
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the photoproduction limit in Q. m = 4.75 (GeV/c 2 ) i.e. the corrected version of 
fig. 11(b) of 

Upon integration of the eqs. ( 3.2S| ), and (3.41) we find agreement with all the 
plots shown in [Q] except for discrepancies with three figures, namely, fig. 9(b), fig. 
11(a), and fig. 11(b). The correct plots are reproduced here as fig. 6, fig. 7(a), and 
fig. 7(b), respectively. This checking also uncovered minor errors in the previous 
FORTRAN programs for the inclusive calculations which have now been corrected. 
Fortunately, they were all in the virtual-photon-quark channels so were not signifi- 
cant numerically and do not alter any of the other plots or results in references QJ, 
01, and 0. 



4.2 Hadronic structure functions and correlations 



Recalling that the probability density is related to the momentum density via 
M/) = A*/) we write (jjj ) as 



(4.5) 



We now specialize to the case where H is a proton, as in the case of HERA. Using 
the relations 

Q 2 



F k 



4vr 2 



(4.6) 



where k = 2,L with 02 = (Jq + 3ai/2, and the relations for the scaling functions 



and (4.4), we find 

Q 2 « s Gu 2 ) 



F k (x,Q z ,m z ) 



4vr 2 



+ 



1 d£ 

Cmin C 

1 di 

£min C 



e 2 H fg/p(C,^ 2 )c 



(0) 



eHfg/p(t,V A ) \ c { kl + C k}g 
,2 



+ E fi/pfo 



i=q,q 



P 2 Jl) , =(1) , n M 



k,i 



(1) 



(4.7) 



where k = 2,L. We have set /if = fj, and shown the dependence of the running 
coupling a s explicitly. The lower boundary on the integration is given by £ m ; n = 
x(4m 2 + Q 2 ) /Q 2 . This formula yields the standard heavy quark structure functions 
F2(x,Q 2 ,m 2 ) and Fl(x,Q 2 ,m 2 ) for electron proton scattering, and we will present 
results as differentials of these functions. 

From the formalism described in the previous section we are left with events 
of positive weight and counter-events of either positive or negative weight. Our 
program outputs the final state four vectors of the event and counter-event together 
with the corresponding weight. We histogram these into bins to produce differential 
distributions. 



14 



We note that the inclusive distributions in the heavy quark transverse momen- 
tum, ( i.e., dF2(x, Q 2 , m 2 , pt) / dpt and dF^^x, Q 2 , m 2 , pt) / dpt) and in the heavy quark 
rapidity ( i.e., dF 2 (x,Q 2 ,m 2 ,y)/dy and dFi(x,Q 2 ,m 2 ,y)/dy) were already pub- 
lished in [J2J for a range of x values at fixed Q 2 and for a range of Q 2 values at fixed 
x. These ranges are covered by the HI and ZEUS detectors in HERA. We have 
reproduced these plots with our new programs. Also, we have previously published 
distributions in the invariant mass of the QQ pair, ( i.e., dF 2 (x,Q 2 ,m 2 ,M)/dM and 
dFL(x,Q 2 ,m 2 ,M)/dM) for charm production in ||. Therefore we concentrate on 
new results involving distributions which are sensitive to the four momenta of all 
the final particles. 

We start by considering the production of charm quarks. We use m = m c = 
1.5GeV/c 2 and simply choose the factorization (renormalization) scale as fi 2 = 
Q 2 + 4(m 2 + (P t c + Pf 5 ) 2 /4). Note that there are many possible choices of scale as 
we have all components of the final four vectors. Aside from the Pt dependence, 
this choice reduces to the usual choice of fi 2 = Q 2 for electroproduction of massless 
quarks and jj, 2 = Am 2 for the photoproduction of charm quarks. We introduce a Pt 
dependence by adding in the average of the magnitude of the transverse momenta 
of the heavy quark and heavy antiquark. As mentioned earlier we use the CTEQ3M 
parton densities in the MS scheme and the two loop a s with A4 = 0.239 GeV. 

The first distribution we present basically measures the transverse momentum of 
the additional jet which recoils against the heavy quark pair. The Pt distribution of 
the charm-anticharm pair is shown in fig. 8(a) where we plot dF 2 (x,Q 2 ,m 2 ,P t )/dP t 
as a function of Pt. The histograms are presented at fixed Q 2 = 12 (GeV/c) 2 for x 
values of 4.2 x 10~ 4 , 8.5 x 1CT 4 , 1.6 x 1CT 3 and 2.7 x 10~ 3 respectively. One sees 
that the Pt distribution peaks at small Pt and has a small negative contribution in 
the lowest bin. This is a region where the dominant contribution is from counter- 
events so the weights can be negative. The results of this calculation require missing 
contributions from even higher order perturbation theory before this bin will have a 
positive weight. In general there will be corresponding bins in all the exclusive plots 
where the weights can be negative. The depth of the negative bins is a function 
of x, Q 2 , and the choice of scale. Note that at larger Pt the structure function 
is dominated by the contribution from the square of the bremsstrahlung graphs 
so the weights are positive. Figure 8(b) shows the corresponding results for fixed 
x = 8.5 x 10- 4 plotted for the Q 2 values of 8.5 (GeV/c) 2 , 12 (GeV/c) 2 , 25 (GeV/c) 2 
and 50 (GeV/c) 2 respectively. The distributions peak near small Pt and are either 
small or negative in the first bin. The histograms with Q 2 = 12(GeV/c) 2 and 
x = 8.5 x 10~ 4 (the dotted line) are the same in figs. 8(a) and 8(b). We have also 
used the same scales on the axes so one can easily see that there is a greater change 
if we fix x and vary Q 2 than if we fix Q 2 and vary x. We will continue to use the 
same scale for all the pairs of later plots to simplify the comparison between them. 

We now turn to the distributions in the azimuthal angle between the outgoing 
charm quark and charm antiquark which we denote as A(p. This is the angle between 
the Pt vectors of the heavy quark-antiquark in the c. m. frame of the virtual- 
photon-hadron system. Since we integrate over the azimuthal angle between the 
plane containing the incoming and outgoing leptons and the plane defined by the 
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incoming parton and outgoing heavy quark (to limit our discussion to F2 and Fl) 
we can only plot relative azimuthal correlations. In the Born approximation this 
distribution is a delta function at tt, as their four momenta must balance. Due to the 
radiation of the additional light mass parton, the distribution has a tail extending 
below tt and has a valley at tt. The distributions become negative in the highest 
bins. This negative region is a general feature of all exclusive calculations. Figure 
9(a) shows results for dF2(x, Q 2 ,m 2 , Acp) / d(Acj)) at the same values of fixed Q 2 and 
variable x as chosen previously in fig. 8(a), while fig. 9(b) shows the results for fixed 
x and variable Q 2 as chosen previously in fig. 8(b). Note again that the dotted 
histograms are the same in fig. 9(a) and 9(b), and there is more variation for fixed 
x and changing Q 2 than for fixed Q 2 and changing x. 

Finally we show distributions in the so-called heavy quark cone size variable 
R = \J (Acj))' 2 + (Ar/) 2 where A<f> is the azimuthal angle between the Pt vectors of 
charm-anticharm quarks in the c. m. frame of the virtual-photon-hadron system 
and A77 is the difference in pseudo-rapidities of the charm quark-antiquark pair. We 
define pseudo-rapidity to be rj = 1/2 In [(1 + cos0)/(l — cos 9)] where 9 is the angle 
the quark (antiquark) makes with the axis defined in the back-to-back photon- 
hadron system. In fig. 10(a) we choose fixed Q 2 and variable x as in fig. 8(a) and 
plot dF2(x, Q 2 ,m 2 , R)/dR. In fig. 10(b) we choose fixed x and variable Q 2 as in fig. 
8(b). The influence of the additional radiation causes these distributions to peak 
below R = tt and develop a dip at values above R = tt. The negative bins are 
unphysical and should be partially filled in by higher order corrections. The depth 
of the negative bins is a function of the bin width. If we use a wider bin, the dips 
will be less pronounced. Again, as in the previous cases, we note there is a greater 
change in the histograms when we fix x and vary Q 2 than when we fix Q 2 and vary 
x. 

The next six figures repeat the previous distributions but for Fl(x,Q 2 ,m 2 ) 
rather than F2(x,Q 2 ,m 2 ). The general features are much the same, but the in- 
tegrated results are smaller for Fl(x,Q 2 ,m 2 ) than for ^(x, Q 2 , m 2 ). Figures 11(a) 
and 11(b) show the distributions in the transverse momentum of the pair. Figures 
12(a) and 12(b) show the distributions in the azimuthal angle between the Pj vec- 
tors of the outgoing charm quark and charm antiquark. Then we present in figs. 
13(a) and 13(b) the histograms of the distributions in R. 

For completeness we now repeat the last twelve distributions taking the heavy 
quark to be the bottom quark with mass mj, = 4.75 GeV/c and the heavy-antiquark 
to be the bottom antiquark with the same mass. The renormalization (factorization 
) scale is chosen to be /x 2 = Q 2 + m 2 + (Pj? + P^) 2 /4. Aside from the Pt dependence, 
this choice reduces to the usual choice of // 2 = Q 2 for electroproduction of massless 
quarks and /j, 2 = m 2 for the photoproduction of bottom quarks. Here we used 
the CTEQ3M ( MS ) distributions and the two loop running coupling with A5 = 
0.158GeV. 

The integrated results are smaller for bottom quarks than for charm quarks 
reflecting a decrease in F2 by a factor of approximately 50, and in Fl by a factor 
of approximately 150. We refer the interested reader to || for tables containing 
F2(x, Q 2 , m 2 ) and Fl(x,Q 2 ,m 2 ) for both charm and bottom quarks in the same x 
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and Q 2 bins presented here. Apart from this decrease there are general features 
which are to be expected. In figs. 14(a), 14(b), 17(a), and 17(b) we see that the 
Pt spectra are harder. The A<f> histograms in figs. 15(a), 15(b), 18(a), and 18(b) 
show dominant back-to-back peaking, which is stronger for bottom quarks than for 
charm quarks as the NLO corrections are correspondingly smaller. Finally in figs. 
16(a), 16(b), 19(a), and 19(b) we see that the R histograms are very Q 2 dependent. 

In addition to checking that all plots presented are indeed independent of the 
parameters p and u we have studied the renormalization (factorization) scale de- 
pendence of the Pt, Atf>, and R distributions presented here and the M distributions 
presented in 0. To aid the discussion we present tables 1-4 containing the aver- 
ages of Pt, Acj), R, and M as functions of the renormalization (factorization) scale 
for charm production. The subscripts 2 and L refer to averages over F2 and Fl 
respectively. As customary we have presented averages for the renormalization (fac- 
torization) scale fj, = fj,Q chosen in the previous plots and for fj, = fJ>o/2 and \x = 2/iq. 
In tables 5-8 we show the same quantities for bottom production. Typical scale 
variations from the central value are around 5 percent for charm production and 1 
percent for bottom production. In general the shapes of the Pt, A(j), and R distri- 
butions remain fixed while we vary the scale but the normalization changes. This is 
to be expected because the corresponding plots are delta functions at lowest order. 
However, the M distributions have contributions from LO already so while the nor- 
malization stays roughly fixed the shape of the plot changes slightly reflecting the 
change in the average values in the tables. 

Our programs can also produce two dimensional plots. For example one might 
study dFk(x, Q 2 ,m 2 ,yQ, y®) / dy® dy® or dFk(x, Q 2 ,m 2 ,yQ, M)/dy®dM or other com- 
binations. In addition, experimental cuts can be implemented. 

5 Conclusion 

In this article we have outlined the NLO calculation of the virtual-photon-parton 
(Wilson) coefficient functions in the exclusive production of heavy quarks plus one 
jet. This completes the study of the NLO electroproduction of heavy quarks. The 
single particle inclusive calculation has already been published [Q], Q. Also there 
are fits available in || for the 77 and £ dependence of the scale independent coefficient 
functions allowing for a fast numerical estimation of integrated rates for experiments. 
In this series of papers we have concentrated on deep-inelastic electroproduction 
where there is no need to introduce any partonic densities in the photon. 

The computer program we have written for the exclusive calculation has the 
advantage that the four vectors of the heavy quark, heavy-antiquark and/or one 
additional light parton jet are produced for each event and can be subjected to 
experimental cuts. We have not done this in any of the plots shown here but the 
computer program is available and can be easily modified to incorporate acceptances 
of the detectors at HERA. [] We have previously presented NLO distributions in the 
invariant masses of the heavy quark heavy-antiquark pair || and shown that they 

1 Requests for the computer program should be sent to smith@elsebeth.physics.sunysb.edu. 
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have reasonably smooth i^-factors, i.e., one can generate these distributions by 
multiplying the Born differential cross sections by constant factors, usually around 
K = 1.3 for bottom and 1.7 for charm. However, the study Q of the single particle 
inclusive distributions in pt and y showed that they do not have smooth K factors. 

For this paper we have presented those plots which depend on information from 
the four vector of the additional jet. We showed the distributions in the transverse- 
momentum (P t ) of the heavy quark antiquark pair, in the azimuthal angle (A0) 
between the Pj vectors of the heavy quark and heavy antiquark, and in the dis- 
tribution in the heavy quark cone size (R). All quantities were predicted in the c. 
m. frame of the photon-proton system after integration over the azimuthal angle 
between the plane containing the incoming and outgoing lepton and the plane con- 
taining the incoming proton and outgoing heavy quark. The results were presented 
as distributions in ^(x, Q 2 ,m 2 ) and Fl(x, Q 2 , m 2 ) at specific points in x, Q 2 and 
m 2 = m 2 or m 2 = m 2 . None of these distributions can be reproduced by any 
-KT-factor multiplication as the corresponding Born distributions are proportional to 
delta- functions. The spread of the tails in these distributions indirectly measures 
the relative size of the NLO contribution to that of the LO contribution. In all cases 
the histograms of these distributions have negative bins. These are regions where 
the NLO calculation is not sufficient and a NNLO order calculation (or some form 
of resummation) should be made. A general statement about the magnitude of the 
NLO contribution compared to the LO one is difficult to make as the size and sign of 
the corrections may vary strongly between different regions of phase space. However, 
we see that all plots have larger Q 2 variation at fixed x as comparied to varying x at 
fixed Q 2 . By varying the renormalization (factorization) scale we observed that the 
distributions presented here changed in normalization but not in shape while the 
invariant mass distribution presented earlier || kept approximately the same nor- 
malization and had a mild shape change. In our study we also calculated averages of 
various quantities and found the typical variation between central and extreme scale 
choices of 5 percent for charm production and 1 percent for bottom production. 
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Table 1 



Range 


(Pth 


<A0> 2 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.73 


1.82 


1.35 


1.45 


2.01 


2.32 


8.5 xlO" 4 


12 


3.02 


2.04 


1.51 


1.42 


1.97 


2.28 


8.5 xlO" 4 


25 


3.72 


2.59 


1.95 


1.39 


1.91 


2.22 


8.5 xlO" 4 


50 


4.53 


3.24 


2.48 


1.40 


1.88 


2.17 


















4.2 xlO" 4 


12 


3.63 


2.45 


1.80 


1.16 


1.79 


2.15 


8.5 xlO" 4 


12 


3.02 


2.04 


1.51 


1.42 


1.97 


2.28 


1.6 xl0~ 3 


12 


2.52 


1.71 


1.28 


1.63 


2.12 


2.39 


2.7 xlO~ 3 


12 


2.14 


1.46 


1.10 


1.79 


2.23 


2.47 



Table 1. Variation of {Pt)i and (A</>)2 for charm production with //q = Q 2 +4(to;? + 
(P t c + P t c ) 2 /4) for various x and Q 2 values. 
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Table 2 



Range 




(M> 2 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.25 


2.70 


2.93 


7.96 


7.58 


7.31 


8.5 xlO" 4 


12 


2.27 


2.70 


2.94 


8.54 


8.14 


7.87 


8.5 xlO" 4 


25 


2.34 


2.76 


3.00 


10.2 


9.77 


9.45 


8.5 xlO" 4 


50 


2.49 


2.87 


3.11 


12.3 


12.0 


11.7 


















4.2 xlO" 4 


12 


2.05 


2.55 


2.83 


8.98 


8.54 


8.21 


8.5 xlO" 4 


12 


2.27 


2.70 


2.94 


8.54 


8.14 


7.87 


1.6 xl0~ 3 


12 


2.45 


2.83 


3.03 


8.09 


7.75 


7.51 


2.7 xlO~ 3 


12 


2.59 


2.92 


3.10 


7.72 


7.41 


7.21 



Table 2. Variation of (R)2 and (M)2 for charm production with ^ = Q 2 + 4(m 2 + 
(P t c + P t c ) 2 /4) for various x and Q 2 values. 
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Table 3 



Range 


(Pth 




X 


Q' 1 


H = no/2 


M = Mo 


H = 2fi 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.47 


1.73 


1.31 


2.17 


2.45 


2.62 


8.5 xlO" 4 


12 


2.78 


1.96 


1.50 


2.14 


2.42 


2.59 


8.5 xlO" 4 


25 


3.58 


2.60 


2.00 


2.12 


2.38 


2.55 


8.5 xlO" 4 


50 


4.60 


3.44 


2.68 


2.13 


2.36 


2.53 


















4.2 xlO" 4 


12 


3.25 


2.30 


1.75 


2.02 


2.33 


2.52 


8.5 xlO" 4 


12 


2.78 


1.96 


1.50 


2.14 


2.42 


2.59 


1.6 xl0~ 3 


12 


2.37 


1.68 


1.28 


2.25 


2.51 


2.66 


2.7 xlO~ 3 


12 


2.05 


1.45 


1.12 


2.34 


2.57 


2.71 



Table 3. Variation of (P t ) l and (A<j>) l for charm production with /j,q = Q 2 +4(m 2 + 
(Pjr + P t c ) 2 /4) for various x and Q 2 values. 
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Table 4 



Range 


(R)l 


(M)l 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.59 


2.78 


2.90 


6.15 


5.97 


5.86 


8.5 xlO" 4 


12 


2.58 


2.77 


2.89 


6.56 


6.36 


6.23 


8.5 xlO" 4 


25 


2.58 


2.77 


2.89 


7.75 


7.53 


7.39 


8.5 xlO" 4 


50 


2.62 


2.79 


2.91 


9.42 


9.17 


9.04 


















4.2 xlO" 4 


12 


2.01 


2.70 


2.84 


6.80 


6.51 


6.34 


8.5 xlO" 4 


12 


2.58 


2.77 


2.89 


6.56 


6.63 


6.23 


1.6 xl0~ 3 


12 


2.65 


2.83 


2.94 


6.39 


6.21 


6.11 


2.7 xlO~ 3 


12 


2.71 


2.87 


2.97 


6.21 


6.09 


5.99 



Table 4. Variation of (R)l and (M)l for charm production with /j,q = Q 2 + 4(m 2 + 
(P t c + P t c ) 2 /4) for various x and Q 2 values. 
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Table 5 



Range 


(Pth 


<A0> 2 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


3.52 


2.54 


1.94 


2.29 


2.52 


2.66 


8.5 xlO" 4 


12 


3.88 


2.83 


2.17 


2.24 


2.48 


2.63 


8.5 xlO" 4 


25 


4.64 


3.47 


2.70 


2.17 


2.41 


2.56 


8.5 xlO" 4 


50 


5.36 


4.11 


3.25 


2.13 


2.35 


2.51 


















4.2 xlO" 4 


12 


4.70 


3.48 


2.68 


2.12 


2.37 


2.54 


8.5 xlO" 4 


12 


3.88 


2.83 


2.17 


2.24 


2.48 


2.63 


1.6 xl0~ 3 


12 


3.16 


2.27 


1.74 


2.35 


2.57 


2.70 


2.7 xlO~ 3 


12 


2.57 


1.82 


1.40 


2.45 


2.65 


2.76 



Table 5. Variation of {Pt)2 and (A<j>)2 for bottom production with = Q 2 + m 2 + 
(P t 6 + P t fe ) 2 /4 for various x and Q 2 values. 
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Table 6 



Range 




(M> 2 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.77 


2.96 


3.07 


16.7 


16.6 


16.4 


8.5 xlO" 4 


12 


2.74 


2.93 


3.05 


17.3 


17.2 


17.0 


8.5 xlO" 4 


25 


2.71 


2.89 


3.02 


18.8 


18.6 


18.4 


8.5 xlO" 4 


50 


2.71 


2.88 


3.01 


20.5 


20.3 


20.1 


















4.2 xlO" 4 


12 


2.64 


2.85 


2.98 


18.3 


18.1 


17.8 


8.5 xlO" 4 


12 


2.74 


2.93 


3.05 


17.3 


17.2 


17.0 


1.6 xl0~ 3 


12 


2.83 


3.00 


3.11 


16.4 


16.3 


16.1 


2.7 xlO~ 3 


12 


2.90 


3.06 


3.15 


15.6 


15.5 


15.4 



Table 6. Variation of (R)2 and (M)2 for bottom production with (Xq = Q 2 + m 2 + 
(P t fe + Pt) 2 /4 for various x and Q 2 values. 
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Table 7 



Range 


(Pth 




X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.89 


2.24 


1.80 


2.65 


2.75 


2.82 


8.5 xlO" 4 


12 


3.29 


2.54 


2.03 


2.62 


2.72 


2.80 


8.5 xlO" 4 


25 


4.08 


3.18 


2.55 


2.56 


2.68 


2.76 


8.5 xlO" 4 


50 


4.86 


3.85 


3.13 


2.52 


2.64 


2.73 


















4.2 xlO" 4 


12 


3.91 


3.05 


2.46 


2.55 


2.67 


2.75 


8.5 xlO" 4 


12 


3.29 


2.54 


2.03 


2.62 


2.72 


2.80 


1.6 xl0~ 3 


12 


2.72 


2.08 


1.67 


2.67 


2.77 


2.84 


2.7 xlO~ 3 


12 


2.25 


1.72 


1.38 


2.72 


2.81 


2.87 



Table 7. Variation of {Pt)L and l for bottom production with = Q 2 + m 2 + 
(P t 6 + P t b ) 2 /4 for various x and Q 2 values. 
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Table 8 



Range 


(R)l 


(M)l 


X 


Q' 1 


H = no/2 


M = Mo 


M = 2 Mo 


M = Mo/2 


M = Mo 


M = 2 Mo 


8.5 xlO" 4 


8.5 


2.91 


2.98 


3.02 


15.0 


14.9 


14.8 


8.5 xlO" 4 


12 


2.88 


2.95 


3.00 


15.3 


15.1 


15.0 


8.5 xlO" 4 


25 


2.84 


2.92 


2.98 


16.0 


15.9 


15.8 


8.5 xlO" 4 


50 


2.82 


2.90 


2.96 


17.1 


17.0 


16.8 


















4.2 xlO" 4 


12 


2.84 


2.91 


2.97 


15.7 


15.5 


15.3 


8.5 xlO" 4 


12 


2.88 


2.95 


3.00 


15.3 


15.1 


15.0 


1.6 xl0~ 3 


12 


2.92 


2.99 


3.04 


14.8 


14.8 


14.7 


2.7 xlO~ 3 


12 


2.96 


3.02 


3.06 


14.4 


14.4 


14.3 



Table 8. Variation of (R) l and (M) l for bottom production with ^ = Q 2 + m b + 
(P t 6 + i- , t b ) 2 /4 for various x and Q 2 values. 
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Appendix A 

Here we discuss the kinematic variables and phase space in the 2^2 and 2^3 
reactions. We then define the generalized plus distributions used in the text. The 
Mandelstam invariants for the reaction, 

1 *( q ) + a 1 (k 1 )^Q{p 1 ) + Q{p 2 ), (A.l) 

with k\ = and p\ = m 2 are 

s' = s - q 2 = (q + fci) 2 -q 2 = 2q ■ k u 

h = t - m 2 = (k\ - p 2 ) 2 - m 2 = -2ki ■ p 2 , 

ui = u - m 2 = (q -P2) 2 - m 2 = -2q ■ p 2 + q 2 , (A. 2) 

which satisfy s' + 1\ + u\ = 0. The two body phase space in space-time dimension 
n = 4 + e is 

where f3 = \/l — p, p = 4m 2 /s and 9\ is the angle between q and p! in the 7*ai 
center-of-mass frame. Therefore we have 

h = --s'(l-/3cos0i), 

U! = -^s'(l + PcosOi), (A.4) 

with < 61 < 7T. 

For the two to three body process 

7*(<?) + ai(fci) Q(Pi) + Q(P2)+a 2 (k 2 ), (A.5) 

with k 2 = and = m 2 there are 5 independent invariants which we take to be 

s' = s -q 2 = (q + A4) 2 - q 2 = 2q • k u 
h = {h -p 2 ) 2 - m 2 = -2ki -p 2 , 
u x = ui-q =(q-p 2 ) -m -q = -2q ■ p 2 , 
t' = fa - k 2 f = -2k x ■ k 2 , 

u > = (q-k 2 ) 2 = q 2 -2q-k 2 , (A.6) 
where q + k\ = k 2 + p\ + p 2 . The remaining five invariants are 

«3 = (k 2 +P2) 2 -m 2 = -s' — u' — t' — u[ —ti = 2k 2 ■ p 2 , 

54 = (k 2 +pi) 2 - m 2 = s' + u\ + ii = 2k 2 ■ p±, 

55 = (Pi + P2) 2 = s' + v! + t' = 2m 2 + 2pi • P2, 
ue = (k\ — pi) 2 — m 2 = —s' — t' — t\ = —2k\ ■ p±, 

u 7 = (q ~ Pi) 2 ~ m2 = ~ s> — u' — u'i = q 2 — 2p\ ■ q . (A. 7) 
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We introduce the variables x = s' 5 /s' where s' 5 = s 5 — q 2 , and y is the cosine of 
the angle between q and k2 in the 7*01 center-of-mass frame. These have ranges 
p* < x < 1 and -1 < y < 1 with p* = (Am 2 - q 2 )/s'. We then find 

t' = -V(^)(l-x)(l + y), 

u' = 9 2 (^)-^(7)(l-x)(l- y ). (A.8) 

In the center-of-mass system of the outgoing heavy quark antiquark pair we decom- 
pose the momenta as follows: 

q = ((A0,0,|q|), 

k\ = 0, sin i/;, cos V>)> 

^2 = (k®, 0) ^1 sin^, I q| +k® cos ip), 

Pi = ^v 7 ^! 1 '^ sin 6*2 sin 6*i, /3 5 cos 6*2 sin 6*i, /3 5 cos 6*i ), 

P2 = ^ v / ^5( 1 '-/ ? 5sin6'2sin6'i,-/35cos6'2sin6'i,-/35cos0i), (A.9) 



where 



|q| 



^2 



s + v! 



2, 


/S5 ' 






2 A 




•55 


-u> 


2 




s • 


" S 5 



cos ip 



-s' + 2k^q° 



2fc°|q| ' 

/3 5 = yjl - Am 2 /s 5 . (A. 10) 

Note that in the limit x — ► 1, — > /?. The remaining two independent invariants 
are 

ii = — -(55 — + /?5 cos 6*2 sin^i sin^ + /?5 cos 6*i cos 
ui = g 2 - i(s + n' + 2^/3 5 |q| cos^i). (A.ll) 
We will also need the following expressions 

s 3 = - [s - s 5 + («5 - u')@5 cos 02 sin 6»i sin^ + (3 5 cos6'i(2 v /iJ | q| +(s 5 - u') cos ip)] , 
S4 = - [s — s 5 — (s 5 — nO^cos^sin^isin^ — /3 5 cos 0i(2y / iJ |q| +(s 5 - u')cosip)] . 
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The three body phase space in n = 4 + e space-time dimensions expressed in 
terms of x,y,6\, and #2 is 



dT 3 = HNdT 



(5) (s>y 



Al+e/2 /o'\ 1 + e / 2 



2tt 



(1 - x) 1+e (l - y 2 ) t / 2 dysm t 9 2 de 2 , (A.12) 



with < 6>i < 7T, < 2 < 7T, p* < x < 1, -1 < y < 1, and 



H 

N 

dT 



r(l + e/2) 



7T 



r(l-e/2)r(l + e) 1 i2 £2 + 0(e3) ' 
(4vr)- e / 2 , . . 1 

-V r(1 - E/2 » = lir,n 2 



(5) 



(4vr) 

22 /85 

16vr Uvr 



e/2 



l+e 



1 



sin 1+e OxdOxdx, 



(A.13) 
(A.14) 

(A.15) 



5 r(l + e/2) 

where 2/e = 2/e + 7^ — ln4-7r. 

The generalized plus distributions encountered in sec. 3 arise by making the 
replacements 



(i-y 2 )~ 1+t - 



1 



l + y 



1 + 2/ 



M 1 +//) (1 + lnw) +(>«•} 



+ 



(A.16) 



1-2/ 



+ [5(1 + y) + 5(1 - y)} ( 1 + ~ ln2u, ) + 0(e) 



(1-x)" 



-l+e 



1 



,1 — X 

+ 5(1 -x 



1 



2e 2 

ln(l — x) 
1 — x 

+ 2 In j3 + 2e In 2 (3 



(A.17) 



+ 0(e 2 ), 



(A.18) 



inside integrations over smooth functions. We have defined $ = y/1 — p, and the 
generalized plus distributions are defined by 



dxf(x) 



1 



1 — x 



dxf(x) 

-1+0) 



/ln(l - x) 
\ 1-x 



1 — X 

1 J _/(*)-/(!) 



<ix ; 



1 - x 



ln(l — x) 



dyf{y) 



l 



dyf(y) 



l + y 

1 



-»y(y)-/(-i), 

1 1 + 2/ 

l-w 1^2/ 



(A.19) 
(A.20) 
(A.21) 
(A.22) 



where p* < /5 < 1 and < cu < 2. Thus we see, in distinction from the phase space 
slicing method, that the soft and collinear parameters p and u are not required 
to be small. It is understood that when the integration range does not enclose a 
singularity, the distribution sign is dropped (c.f. ([jug)). 
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Appendix B 

In this appendix we discuss the derivation of the soft limit x = 1 of the matrix 
element Mf (3). We use standard techniques for the emission of soft gluons [28], 

Consider the diagram shown in fig. 3(d). The momenta are defined according 
to fig. 1. Take the incoming gluon to have color index a and Lorentz index /i, the 
outgoing gluon to have color index c and Lorentz index a, and the photon to have 
Lorentz index p. Using standard Feynman rules one finds 



Lk\ ■ K 2 



(B.l) 



where 



Vfxraiki, -k 2 , k 2 - fei) = (h + k 2 ) a g^ a + (h - 2k 2 ) fM g CTa + (k 2 - 2k{] a g m . (B.2) 

By is the sum of the lowest order diagrams in fig. 2 and we have chosen the Feyn- 
man gauge. Now by using current conservation 

and the transversality of the gluon polarization tensor [e^(ki)ki = 0], and recalling 
that in the soft limit k 2 — > 0, a short computation shows that 



k 2 ) a B ba P(q,k 1 -k 2 )=0 



igf 



abc 



fci-ifea " 



e At (/ci) J B bw (g, fci) + finite terms as k 2 -> 0. 



(B.3) 



As we work in the Feynman gauge we must consider the interference between all 
diagrams in fig. 3 |32|. Analyzing them as we have fig. 3(d) leads to 

Mf (3) = %g A e 2 e 2 H NC F [C A S i>OK + 2C F Si, QED ] (B.4) 

in the soft limit when one neglects terms that are finite when multiplied by the 
(1 — x) 1+6 term in the three body phase space. We have defined 



Si,OK 
Si,QED 



[(kip 2 ) + (pih) - (pip 2 )] B it QED 

{pm) - tKpipi) 



SP2P2) 



B 



i,QED- 



(B.5) 



and denoted (uv) = u ■ v/(u ■ k 2 v ■ k 2 ) which is often called the eikonal factor. 
Expanding the factors (uv) we reproduce eqs. (3.21) and (3.22) of Recalling the 
definition of ff{0\) we find 



f!(0i) 



f\2 



8g 4 e 2 e 2 H NC F 



B 



i,QED 



CA I (hups) + hviki) ~ ^(piP2) + 2< ^F 



(P1P2) 



2 (PlPl) 



2^(P2P2) 



(B.6) 



with 



dy / sin e 6 2 d6 2 (uv){l - x) 2 (l - y 2 ) 



x=l 



(B.7) 
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Using the expansion of S3 and S4 about x = 1 



S3 

Si 



-s'(l-x) 
2 v ; 



1 + (5 J 1 — y 2 cos 62 sin 9\ + fly cos 9\ 



\s'(l-x) 



1 — (3y 1 — y 2 cos 6*2 sin 9\ — fly cos 8\ 
we can write the integrals I(u V ) in terms of 



+ 0((l-xf), 
+ 0((l-x) 2 ), 



1^=1 da(sma) n - 3 dfl(sin fl) n ~ 4 (A + B cosa + C sinacos fl)~ j , (B.8) 



with j = 1,2 and A 2 / B 2 + C 2 , and 



7 



da(sina)"- 3 / d/3(sin/3) n - 4 (A+S cos a+C sin a cos fl) 3 ' (a+b cos a) 



(B.9) 

-6, and i = j = 1. The second integral may be found in 
the literature [29 and the first evaluated using the methods of [3C]. We find 

'a + vb 2 + c 2 \ 



with A 2 ^ B 2 + C 2 , a 
liter: 



7T 



In 



VB 2 + C 2 j \A-VB 2 + C 2 
-(n-4) 



Li 2 



2V^ 2 + C 2 \ 1 2 /i + VF+?J2 
' - In 1 



A + y/B 2 + C 2 / 4 [a-VWTc 2 



(B.10) 



and 



r(2) 



2tt 



^2 _ ^2 _ C 2 



1, A A + VB 2 + C 2 

1 (n - 4) ; = In ; 

2 K ' ^B 2 + C 2 \ A - VB 2 + C 2 



(B.ll) 



where we drop 0((n — 4) 2 ) terms and quote, for completeness, the result for the 
second integral M 



(1,1) 



7T 



2 / (A + B) 2 

+ In ' 



a(A + B)\n-A \A 2 - B 2 - C 2 



hi 2 



A- VB 2 + C 2 \ l-fA + ^WTC 2 

In 

/ 2 



+ 2Li 2 



A + B 

b + Vb 2 Tc 2 \ 



A - VB 2 + C 2 



2Li, 



'g- VB 2 + C 2S 
A + B 



(B.12) 



again dropping 0((n — 4) 2 ) terms. The dilogarithmic function Ui2(x) is defined in 
Prom these integrals and (PH), (|B~5p, (^6|), and (U) w e find 



dag = C i>g M* oft dr 2 



(B.13) 
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with 



M, 



soft 



„4 2 2 



e e H N/j, e C e CF CaS^ok + ^CfSi^qed 



B 



i,QED- 



(B.14) 



The factor 



1 / 2\- e / 2 

^_ p e(7 B -ln47r)/2 ^_ 



(B.15) 



is common not only to da\ s J but to da^g as well (see Appendix A of |p ). The 
mass parameter fx originates from the dimensionality of the gauge coupling in n 
dimensions. The remaining terms are 



S, 



OK 



e 2 e 



1 / 2m 2 



1 In 



1+0 



41n/3- In— ^ - In 



-Til 



In At + 2 In + 2 In — ^ - 4 In 

m z m z m z m z 



In 



In — T + In — ^ In — ^ 

m z m z m z 



In 



1 / 2m 2 



m 2 f3 



81n/3- 21n^ + 41n^ + In 



1-P 



m 2 m 2 ' 1 + (3 



+ > ( In + In — ^ 



ln rrl +4Li2 

ln/3 + 16 In 2 /3 - 3((2) - In 2 



2/3 
1 + /3 



+ 2Li 2 1 + 



2ti 



2ui 



2Li 2 (l + — (1 + /3) 



+ 2Li2 i 1 + 7(1-/5) 



2Li 2 (l + — (1 + /3) 



+ In 2 



In 2 



(B.16) 



and 

Si,QED 



1 / 2m 2 



1 In 



1-/3 
1 + 13 



+ 



1 /2m 2 



/3 



8 In $ - 2 In —= + 4 In + In 



77 V - 



81n/3 + 21n — -41n — - - In 



1-/3 
1 + /? 
1-/3 



In 



1-/3 
1 + /? 



+ 4Li 2 



2/3 
1 + /3 



m- 



m 2 /3 1 + /3' 



(B.17) 



where C(2) = tt 2 /6. 
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Appendix C 



In this appendix we discuss the derivation of the collinear limit y = — 1 of the 
matrix element Mf(3). A similar analysis follows for Mf(3) for which we simply 
quote the results at the end. We closely follow appendix B of thereby deriving 
ff{x, —1,61,62) in n dimensions (see also appendix C of (jOJ). 

Consider the diagrams shown in fig. 3. The momenta are defined according to 
fig. 1. Take the incoming gluon to have color index a and Lorentz index fj,, the 
outgoing gluon to have color index c and Lorentz index a, and the photon to have 
Lorentz index p. Then the sum of the diagrams is 



M ac P = e ^ki)e tr (k 2 ) 



(C.l) 



where D^f is the contribution from fig. 3(d) and BfjEf represents the remaining 
terms that are regular when y = — 1. Using standard Feynman rules one finds 



r> ac P 

fJUT 



with 



q,k 1 -k 2 ), (C.2) 



V flaa (ki,-k 2 , k 2 - fa) = (fa + kfiagftcr + (fa - 2A; 2 ) M 5 CTQ + (k 2 - 2fa) a g fia , (C.3) 

and B baf} the sum of the lowest order diagrams in fig. 2. We choose the propagator 
in the light-like axial gauge where 



P^u(k) = -g^ v + 



kpj]u + k u rj^ 
77 • k 



(C.4) 



Now r] and ki define a plane. Take k± perpendicular to this plane and choose r\ 
such that 7] ■ k\ 7^ and rj 2 = 0. Then decomposing k 2 = (1 — x)fa + r]^ + k± 
implies that £ = —fa\/[(l — x)2r\ ■ fa] and t' = kj_/(l — x). The collinear limit 
we require is k± — > with x fixed. If we use e(ki) ■ fa = 0, e(k 2 ) ■ k 2 = 0, and 
(fci - k 2 ) a P$(fa - k 2 ) = 0(k\) then 



V" a (ku -fa, fa - fa) ~ -klsT - 2^5™ + z + ■■■ + o(k 2 ± ) , (C.5) 



X 



where • • • represents terms that vanish upon contraction with the polarization ten- 
sors. Hence ( |C.1| ) becomes 



M acp 



-igf 



acb 



2 2 , 2 
-k± a gp a ~ 2k±^g aa + - - - k± a g^ a H h 0(k A 



-5 bd P%(fa 



(fa - fc 2 ) 2 



5<^(g, fa - k 2 ) + Jtfg? e»(fa)?(k 2 



(C.6) 



From this result one can see that M ~ 0(1/ yt' ) + R + O(t') as i' — > 0. Hence when 
we square M we can drop interference terms containing R because only the square 



33 



of the first term gives a 1/t' singularity. Now because ki — = xk\ + 0(k±) we 
have B al3p (q, k\ — k2) = B a @ p (q, xk\) + 0(k±) so we can write 



M acp ( M 



acX 



^^B a ^(q, xh ) (B aaX (q, xh )) f P af3 (h 



(l-x) + 



1 1 
+ - 



+ 



^/r^,.,-/,- l ;(/r'\,.,-/, r .^ i 

k± a k±p 



A: i 



+ 



P. 



a/3 



n 



(C.7) 



Using the definition of ff(x, y, Q\, 62) in eq. (3.3) one can easily make the identifica- 
tion used in the main text, namely that 



fi(x, -1, 9 1 ,9 2 ) = ff(x, 6,) + ft (x, 0i, 2 ). 
The second term vanishes upon integration over 02 and 

(*0 2 



(C.8) 



5127rV- 2e a2 e 2 e^iVCA<V 



1 



x(l — x) + 



1 



+ 



1 — x 



Bi,QED (xk\). (C.9) 



The mass parameter \jl originates from the dimensionality of the gauge coupling in 
n dimensions. For the quark channel a similar analysis holds with the result 



ff{x, -1, 0i, 2 ) = f!(x, 0i) + f9(x, 0i, 2 ) 
where again the second term vanishes upon integration over 02 and 



(CIO) 



-l2M 2 ^- 2€ a 2 s e 2 e 2 H NC F {\ + e/2)" 1 
1 + (1 - x) 2 + ex 2 /2~ 



B. 



i,QED 



(xki). 



(C.ll) 



Note that we need the full ff(x, — 1, 0\, 62) in n = 4 dimensions in the plus 
distributions used in the main text not just ff{x,Q\). For this one may take the 
n = 4 limit of eq. flCij ) by using an explicit representation of k±. However, we 
found it instructive to take the collinear limit of Mf(3) in the Feynman gauge and 
compare with the axial gauge results. It is well known [|32| that if one chooses a gauge 
different from the axial one in QCD, then interference graphs also become singular 
in the collinear limit. The axial gauge result follows from the previous analysis. In 
the Feynman gauge the squared matrix element arising from the diagrams in fig. 3 
(including interference terms) has the form 



M '( 3 )=(^2+f + C + 0(t'), 



(C.12) 
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where A, B, and C are functions of the other nine Mandelstam invariants defined 
in appendix A. Expanding the invariants u\, uj, and t\ about t' = we find 



Ui 

h 



Xt\ -j 



«i + Q 2 /?5 cos 6*1 



x^ - Q 2 /? 5 cos Q\ 



t\ - —V2at'(3 5 cos 02 sin 6>i 



where 



and 



+ —<Ui + j3 5 cos 6>i 



-2(1 - x) 



x — 1 



•S'5 



(xs') 



^2 



(C.13) 
(C.14) 



-s'(l - ^ 5 cos 6>l) , 
-4(l + /? 5 cos0i). 



(C.15) 



Direct substitution of the invariants ( p.!3| ) into ( C,12| ) and using the definition of 
Jq{x, y, 0i, 62} in eq. (3.3) yields the results 



f G {x-l,0 x ,0 2 ) 



A 2 



+ 



32(2m 2 - Q 2 )(xs') 
{xt\)u\s 
32(x 2 -x + l) 2 (xs') 2 

X 4 S 



g(x,0 1 ,0 2 ) 

' (xs') 2 
{xt\)ui 



f 9 L (x,-l,0x,0 2 



64Q ; 



g(x,0i,i 



(C.16) 
(C.17) 



where 



'i,t/2j 



2(x 2 - 2x + 2) 



x- 



xs 



A2 



(xtf )uj 



•55 



(1 -xf 



S5P 2 cos 2 02 sin 2 6*i 



(xs 



A2 



(C.18) 



Both /§(x, — 1, 0i, 02) and f 9 L {x,—l,0\,0 2 ) have a common factor g i e 2 e 2 H NCACF 
which we have not shown. We have used the relation xs' + xt\ + u\ = freely 
in deriving ( |C. 16| ) and ( C.17| ). A similar analysis for the quark channel gives the 
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results 



(xt\)u\ 
8(x 2 -2x + 2) 



x* 



l\2 



(xs') 
(xt\)u\ 



fl( x ,-i,e 1 ,e 2 ) 



where 



h(x, i 



A2 



lj "21 



m 



{xs') 
(xt\)u\ 



i 1 - X ) o2 2/i -2/i i XS ') 

- — 2 — S50§ cos 62 sin 6\ 



l\2 



(xti)ui 



(C.19) 
(C.20) 

(C.21) 



Both /q(x, — 1, 61,62) and f[(x, — 1, 61,62) have a common factor g 4 e 2 e 2 H NCF which 
we have not shown. We have checked the formulae for the / 4 ai (x, —1,61,62) (cti = 
g,q, i = G,L) presented above by comparing them numerically with the y — ► — 1 
limit of the complete expression / 4 ai (x, y, 6\, 62). 
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Figure Captions 



fig.l. Illustration of the basic reaction ( |2.lD for heavy-flavor production in virtual- 
photon-parton collisions. 

fig. 2. Lowest order Feynman diagram contributing to the amplitude for the gluon 



fusion reaction (2.12). Another diagram is obtained by reversing the arrows 



on the heavy quark lines. 

fig. 3. The order eg 2 diagrams contributing to the amplitude for the gluon-bremsstrahlung 
reaction ( |3.1| ). Additional graphs are obtained by reversing the arrows on the 
heavy quark lines. 

fig. 4. One loop order eg 3 diagrams for the virtual corrections contributing to the 
reaction (§11). Except for the non-planar diagram additional diagrams are 
obtained by reversing the arrows on the heavy quark lines. The solid line in 
the gluon self-energy graph represents gluon, ghost and quark loops. 

fig. 5. The order eg 2 diagrams contributing to the amplitude for the reaction Q3.30Q . 
Additional graphs are obtained by reversing the arrows on the light-quark lines 
(dashed). 

fig. 6. The r] dependence of the scaling function Cj} q ij]^) in the MS scheme for 
several values of Q 2 (in units of (GeV/c) 2 ) with m = 4.75 (GeV/c 2 ), i. e. 
the corrected version of fig. 9(b) of ref. Q. The solid line corresponds to 
Q 2 = 0.01, the dotted line to Q 2 = 1, the short dashed line to Q 2 = 10, the 
long dashed line to Q 2 = 100, and the dot-dashed line to Q 2 = 1000. 

fig. 7. (a) The rj dependence of the scaling function d^ q (r], £) in the MS scheme for 
several values of Q 2 (in units of (GeV/c) 2 ) with m = 4.75 (GeV/c 2 ), i.e. the 
corrected version of fig. 11(a) of ref. [0. The short dashed line corresponds 
to Q 2 = 10, the long dashed line to Q 2 = 100, and the dot-dashed line to 
Q 2 = 1000. (b) The rj dependence of the scaling function d^ q (r),£) in the MS 
scheme for the same values of Q 2 with m = 4.75 (GeV/c 2 ) i.e. the corrected 
version of fig. 11(b) of ref. [Q]. 

fig. 8. (a) The distributions dF2(x,Q 2 ,m 2 , Pt)/dPt for charm-anticharm pair pro- 
duction at fixed Q 2 = 12 (GeV/c) 2 with x = 4.2 x 10~ 4 (solid line), 8.5 x 10" 4 
(dotted line), 1.6 x 10 -3 (short dashed line) and 2.7 x 10~ 3 (long dashed line), 
(b) The same distributions at fixed x = 8.5 x 10 -4 and Q 2 = 8.5 (solid line), 
12 (dotted line), 25 (short dashed line), 50 (long dashed line) all in units of 
(GeV/c) 2 . 

fig. 9. (a) The distributions dF2(x,Q 2 ,m 2 ,A(j))/d(A.(j)) for charm-anticharm pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distri- 
butions at the x and Q 2 values given in fig. 8(b). 
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fig. 10. (a) The distributions dF2(x,Q 2 ,m 2 ,,R)/dR for charm-anticharm pair pro- 
duction at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 11. (a) The distributions cIFl(x, Q 2 , m 2 , Pt) / dPt for charm-anticharm pair pro- 
duction at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 12. (a) The distributions dF^x, Q 2 , m 2 , A<ft) / d(A<fi) for charm-anticharm pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 13. (a) The distributions dF^x^Q 2 } m 2 , R) / dR for charm-anticharm pair pro- 
duction at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 14. (a) The distributions dF 2 (x, Q 2 , ml, P t )/dP t for bottom-antibottom pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distri- 
butions at the x and Q 2 values given in fig. 8(b). 

fig. 15. (a) The distributions dF 2 (x, Q 2 ,m 2 , Atfi) / 'd(Aip) for bottom-antibottom pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 16. (a) The distributions dF 2 (x, Q 2 ,m 2 , R)/dR for bottom-antibottom pair pro- 
duction at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 17. (a) The distributions dFL(x,Q 2 ,m 2 ,Pt)/dPt for bottom-antibottom pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distri- 
butions at the x and Q 2 values given in fig. 8(b). 

fig. 18. (a) The distributions dFi(x, Q 2 ,m 2 , Acf)) / d(A(f>) for bottom-antibottom pair 
production at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 

fig. 19. (a) The distributions dFi(x, Q 2 ,m 2 , R)/dR for bottom-antibottom pair pro- 
duction at the x and Q 2 values given in fig. 8(a). (b) The same distributions 
at the x and Q 2 values given in fig. 8(b). 
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